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The astrophysical factor for Be(p, 7) B at zero energy, 017(0), is determined from an analysis of Pb( B, 
p± 7 Be) 208 Pb at 52 MeV/nucleon by means of the method of continuum-discretized coupled-channels (CDCC) 
taking account of all nuclear and Coulomb breakup processes. The asymptotic normalization coefficient (ANC) 
method is used to extract Sn(0) from the calculated breakup cross section. The main result of the present paper 
is Sir (0) = 20.9ti'g e ^ b. The error consists of 8.4% experimental systematic error and the error due to the 
ambiguity in the s-wave p- 7 Be scattering length. This value of <Si7(0) differs from the one extracted with the 
first-order perturbation theory including Coulomb breakup by dipole transitions: 18.9 ± 1.8 eV b. It turns out 
that the difference is due to the inclusion of the nuclear and Coulomb-quadrupole transitions and multi-step 
processes of all-order in the present work. The p- 7 Be interaction potential used in the CDCC calculation is also 
used in the ANC analysis of 7 Be(p, 7 ) 8 B. The value of 5 17 (0) = 21.7± ';g eV b obtained is consistent with 
the previous one obtained from a precise measurement of the p-capture reaction cross section extrapolated to 
zero incident energy, 517(0) = 22.1 ± 0.6 (expt) ±0.6 (theo) eV b, where (theo) stands for the error in the 
extrapolation. Thus, the agreement between the values of 517(0) obtained from direct 7 Be(p, 7) 8 B and indirect 
8 B-breakup measurements is significantly improved. 

PACS numbers: 24.10.Eq, 25.60.Gc, 25.70.De, 26.65.-rt 



I. INTRODUCTION 

The solar neutrino problem is one of the central issues in the 
neutrino physics 1 1]. The current interpretation of the solar- 
neutrino deficit is the neutrino oscillation induced by the mass 
difference among v e , v^, and v T , and their mixing angles 
Attention is now focused on the determination of the oscilla- 
tion parameters. The cross section of the p-capture reaction 
7 Be(p, 7) 8 B at incident energy in the center-of-mass (cm.) 
frame of the p- 7 Be system £7 c . m . ~ plays an essential role 
in the solar-neutrino phenomenology, since the observed flux 
of 8 B neutrino is proportional to it; the magnitude of the cross 
section <J pl (E c , m ) is customarily expressed by the astrophys- 
ical factor S\r{E c .m.) = & pl (E c . m .)E c . m _ exp[27r7y], where 
i] is the Sommerfeld parameter. The required accuracy of 
Si7(0), to determine the neutrino oscillation parameters with 
sufficient accuracy, is the error within 5% |3]. 

Recently, precise measurement of cr p7 (i? c . m .) was carried 
out by Junghans et al. 01 at energies of E c , m . =116-2460 
keV, which are low but still higher than stellar energies (~ 20 
keV). Extrapolating the measured Su(E c . m .) to E c . m . = 
using a three-cluster model |5] for 8 B structure, they derived 
517(0) = 22.1 ± 0.6 (expt) ±0.6 (theo) eV b. The three-body 
model, however, did not simultaneously reproduce the mag- 
nitude and the energy-dependence of Sn(E c .m.) sufficiently 
well. Moreover, as pointed out in Ref. @], the uncertainty 
of the s-wave p- 7 Be scattering length (with about 50% error) 
prevents one from determining 5i 7 (0) with very high accu- 
racy. 

Because of the difficulty of the direct measurement of 
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(Tp^(E c , m .) at stellar energies, alternative indirect measure- 
ments were proposed. Coulomb dissociation Q 0, [k], El 
of 8 B is one of such indirect measurements. So far, 
extraction of £17(0) from these experiments has been based 
on the virtual photon theory with the assumption of 8 B dis- 
sociation by virtual electric dipole (El) photon absorption. 
Nuclear interaction and absorption of quadrupole (E2)- and 
multi-photons were not taken into account. The value thus ex- 
tracted from the RIKEN experiment at 52 MeV/nucleon was 
5 17 (0) = 18.9 ± 1.8 eV b H. In the analysis of the MSU 
experiment measured at 44 and 83 MeV/nucleon, 

the E2 contribution to one-step process was estimated from 
the parallel-momentum-distribution of 7 Be fragment and sub- 
tracted from the breakup spectrum of 8 B. As a result, the ex- 
tracted £17(0) was 17.8jli'2 e V b, which is smaller than the 
value obtained at RIKEN mentioned above. However, the 
analysis of the angular distribution of the 8 B breakup cross 
section of the RIKEN experiment showed no contribution of 
E2 transitions. Moreover, the recent experiment at rather high 
energy, 250 MeV/nucleon, done at GSI 1 1 1 ] showed that the 
E2 contribution was negligibly small; the resulting Si 7 (0) was 
18.6 ± 1.2 ± 1.0 eV b. Even though the significance of E2 
transitions can be energy-dependent, the conclusions from the 
MSU and RIKEN measurements at similar energies look in- 
consistent and roles of the E2 component are still not clear. 
More seriously, there exists a non-negligible discrepancy of 
about 15% between the values of Si7(0) mentioned above 
that are derived from direct p-capture and indirect 8 B-breakup 
measurements. 

Very recently, it was shown by a semiclassical calculation 
of 208 Pb( 8 B, p± 7 Be) 208 Pb at 52 MeV/nucleon that the dis- 
crepancy mentioned above was significantly reduced by tak- 
ing account of nuclear-breakup components, E2, and higher- 
order Coulomb breakup processes 1 14]. Motivated by this re- 
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suit, we attempt in the present paper to extract a reliable value 
of 6*17(0) analyzing the 8 B dissociation experiment measured 
at RIKEN 0,|aEi, with the method of continuum-discretized 
coupled-channels fl5l lliill (CDCC), assuming a p+ 7 Be+A 
three-body model of the system, where A stands for the tar- 
get nucleus. The result of the analysis is then used to extract 
5*17(0) by means of the asymptotic normalization coefficient 
(ANC) method |17]. CDCC has been successful in describ- 
ing various processes in which effects of proje ctile-breakup 
are essential OS d S3 ^ E3 |H |M |3 ^ . It has been 
successfull y applie d also to 8 B nuclear and Coulomb breakup 
processes 1 1 3l 1271 1281 12911 . In the CDCC calculation of the 
present work, we include all nuclear and Coulomb breakup 
processes, and take account of the intrinsic spins of both p 
and 7 Be using the channel spin representation and the con- 
sistency of the p- 7 Be interaction potential used in the CDCC 
calculation with the s-wave p- 7 Be scattering length. The wave 
function of 8 B used in the CDCC calculation is found to be 
consistent with the measured cross section of 7 Be(p, 7) 8 B at 
low energies. We then use the calculated breakup cross sec- 
tion and the p- 7 Be separation energy in the 8 B nucleus to ob- 
tain 5i7(0) by the ANC method. An important advantage of 
the ANC method is that there is no restriction of the reac- 
tion mechanisms. In addition, the uncertainty of 5i7(0) due 
to the use of the ANC method can quantitatively be evalu- 
ated 0H. 

The construction of the present paper is as follows. First, 
we describe how to extract the ANC from the three-body 
model CDCC analysis of 8 B breakup. We show that under 
a condition, which is satisfied in the present analysis, the cal- 
culated cross section is proportional to the square of the ANC 
for the ground state of 8 B. This is by no means trivial since the 
CDCC wave function contains not only the ground state of 8 B 
but also its continuum states. We describe in the later section 
the condition and the angular region of 8 B breakup in which 
it is satisfied. Second, we determine 5i7(0) from 208 Pb( 8 B, 
p+ 7 Be) 208 Pb at 52 MeV/nucleon by the ANC method, taking 
account of all order nuclear and Coulomb breakup processes. 
It should be noted that in our CDCC calculation interference 
between nuclear- and Coulomb-breakup amplitudes is explic- 
itly included. We make use of the eikonal-CDCC method 13011 
(E-CDCC), with some modifications, to obtain the scatter- 
ing amplitudes corresponding to large values of the orbital- 
angular-momentum L of relative motion between the projec- 
tile and the target nucleus. We then construct a hybrid ampli- 
tude from the results of E-CDCC and CDCC. This makes it 
possible to carry out very accurate analyses with high com- 
putational speed 1 30]. Uncertainties of the extracted 5i7(0) 
are carefully evaluated. Finally, we discuss the contributions 
of nuclear breakup, E2 transitions, and higher-order processes 
on the extracted value of 5i7(0). The ANC method is also ap- 
plied to direct 7 Be(p, 7) 8 B measurement assuming the same 
p- 7 Be potential in the calculation of the direct capture process 
as used in the CDCC analysis. We show the difference be- 
tween the values of 5i7(0) derived from direct and indirect 
measurements is essentially removed by the present analysis. 

In Sec. Ill Al we describe how to extract the ANC from the 
CDCC analysis of 8 B breakup. Formulation of modified E- 



CDCC and construction of the hybrid scattering amplitude are 
described in Sec. Ill Bl We analyze in Sec. llII Bl the 8 B breakup 
cross section measured at RIKEN with CDCC and extract 
5i7(0) by the ANC method. Uncertainties of the extracted 
5i7(0) are quantitatively evaluated in Sec. llII Cl In Sec. llII PI 
the value of the extracted 5i7(0) is compared with the result 
of the virtual photon theory and the roles of nuclear breakup, 
E2 transitions, and higher-order processes are discussed. We 
discuss in Sec. IIII El the application of the ANC method to 
precise direct measurement of 7 Be(p, 7) 8 B that confirms the 
consistency of its result with the 8 B breakup analysis. Finally 
summary and conclusions are given in Sec. llVI 

II. FORMULATION 

A. Extraction of ANC from 8 B breakup 

In this subsection we describe our method of extract- 
ing 5i7(0) by p+ 7 Be+A three-body model CDCC analy- 
sis of elastic 8 B breakup reaction combined with the ANC 
method 11711 . We recapitulate the ANC method below, partly 
to introduce notations. At very low incident energies, the p- 
capture reaction 7 Be{p, 7) 8 B leading to the ground state of 8 B 
is extremely peripheral because of the Coulomb barrier. The 
T-matrix element of the reaction, therefore, depends on the 
overlap between the ground state wave functions of 7 Be and 
8 B only in the tail region, where the radial part of the overlap 
has the following form 

V4°)(r) = (4 0) (£ 7 )l4 0) (£7,r)) C7 ,r 

~ CW_ vA>+1/2 (2k r)/r, r > R N , (1) 

where i?N is the range of the nuclear interaction between p 
and 7 Be of about 4 fm, (j)^ (4>^) is the ground state wave 
function of 7 Be ( 8 B), k$ is the relative momentum in the unit 
of h between p and 7 Be, £7 stands for the internal coordinates 
of the seven-nucleon system, and W_ v ,£ +i/2 is the Whittaker 
function where £0 is the orbital angular momentum of the p- 
7 Be relative motion in the ground state of 8 B; £q = 1 in the 
present case. The constant C in Eq. Q is the ANC of our 
interest. The cross section <r P7 (0), and consequently 5i7(0), 
are proportional to C 2 . Thus, one can evaluate them once the 
ANC, C, is known. The cross section of any reaction in which 
a proton is transferred to or from 8 B is proportional to C 2 in 
DWBA if the reaction is peripheral in the sense that it takes 
place when the transferred proton is distant from 7 Be more 
than the range i?N- Based on this idea, DWBA analyses of 
proton transfer reactions such as 10 Be( 7 Be, 8 B) 9 Be 13111 and 
14 N( 7 Be, 8 B) 13 C 1 32] have successfully been made to obtain 
the required ANC. 

In the present paper we analyze the data of 8 B elastic 
breakup in which the fragments 7 Be and p are ejected into 
very forward angles. For brevity, we henceforth call such a 
process breakup into forward angles. Since the breakup in 
such a case takes place at large distance R of 8 B from the 
target nucleus, it is mostly due to Coulomb interaction. The 
operator of the interaction is proportional to r™ with n > 1 . 
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Therefore, the dominant contribution to the T matrix is from 
the region of r larger than the range of nuclear interaction 
between the fragments. This is numerically confirmed as de- 
scribed in Sec. IIIICl The reaction, therefore, is peripheral in 
the sense described above. Thus, the ANC can be obtained 
if its cross section is reproduced and is shown to be propor- 
tional to C 2 by theory. We show in the following that the 
cross section is well reproduced by CDCC and is in good ap- 
proximation proportional to C 2 in the present case. 

We describe 8 B breakup at 52 MeV/nucleon with the 
p+ 7 Be+A three-body model illustrated in Fig. 1. Let us as- 




FIG. 1: Illustration of the p+ 7 Be+A three-body system, where A is a 
target. 



sume for simplicity that A is infinitely heavy. The triple dif- 
ferential cross section for ( 8 B, 7 Be+p) reaction is given by 
Vp\1\ 2 with a given constant T> and the three-body phase 
space factor p 13311 . The T-matrix element of 8 B breakup is 



with 



(XiX7<4V|*>, 



V = v A7 (R A7 ,&) + C/ai(Rai) +v 17 (t,&), 



(2) 



where vxy (Uxy) is the interaction (distorting potential) be- 
tween two nuclei X and Y, where the subscripts 1, 7, and A 
stand forp, 7 Be, and A, respectively. The function \i (X7) is 
the plane wave of relative motion between A and the outgoing 
p ( 7 Be), and \& is the exact wave function of the system with 
the incident plane wave of relative motion between 8 B and A 
in the initial channel. 

In the RIKEN experiment that we analyze in the present 
work, emitted fragments were measured at forward angles 
with the 7 Be in the ground state |H 0]. We assume that 
7 Be stays in its ground state in the course of 8 B breakup. 
We have confirmed the validity of this assumption by a 
3 He+ 4 He+ 208 Pb three-body model CDCC calculation of the 
scattering of 7 Be on 208 Pb at 52 MeV/nucleon. The cross sec- 
tions of non-elastic processes including breakup of 7 Be are 
less than about 1 % of that of elastic scattering for the scatter- 
ing angle of the cm. of 7 Be, 87, below about 4°. In general 
67 can differ from the 8 B scattering angle 9 S in the 8 B+ 208 Pb 
— » p+ 7 Be + 208 Pb process, but the difference between the 
two angles is very small, since the energy transfer in the 8 B 



breakup concerned is much smaller than the incident energy 
of 8 B. These facts imply that in T one can neglect matrix el- 
ements of VA7 between the ground state and excited states 
of 7 Be. One can, therefore, approximate va7 by the single- 
folding potential 



C/ A7 (R A7 ) = (4 0) (e7)|«A7(RA7,e 7 )l4 0) (e7)) 57 . 

Similarly, vn can be approximated by 

y 17 (r)EE(4°)(^)|x; 17 (r,6)l4 0) (e7)) e7 . 
In this approximation, T is given by 

I» <XlX7|C/A7 + t/Al 

with 



Vn\0) 



O 



= 1A0) 



(3) 



(4) 



If the incident momentum of the 8 B projectile is P, one can 
make use of the well-known formula for ^ and write 



O 



is 



(5) 



E- H + ie 
where H is the exact Hamiltonian of the system 

H = T r +T R + V + h 7 , 



where T r and Tr are the kinetic energy operators with re- 
spect to r and R, respectively, and h 7 is the internal Hamil- 
tonian of 7 Be. The total energy of the system is denoted by 
E, which is related to the cm. incident energy of 8 B E in by 
E = E- m + e 7 + ei7, where e 7 is the eigenenergy of 4> 7 °^ and 
en is the relative energy between p and 7 Be in the ground 
state of 8 B. Under the assumption of the inert 7 Be nucleus 
mentioned above, one can approximate H in Eq. (|5j by 



H « e 7 + i?3-body, 



(6) 



where -ff3-body is the p+ 7 Be+A three-body Hamiltonian de- 
fined by 

^3-body =T r +T fl + ^ 17 (r) + [/ A 7(RA7) + t/Al(RAl). (7) 

Consequently, 



where 



Uai + U A 1 + Vl7 



n(+y o >(r) e iP - R V (8) 



IS 



E — &J — i?3-body + iS 



(9) 



is Miller's wave matrix with the outgoing boundary condition 
and 



^°)(r) = (^(e 7 )|^(e 7 ,r))e 



,(0) 



(0), 



(10) 



is the overlap function between the ground state wave func- 
tions of 7 Be and 8 B. As one sees from Eq. (|8j, ip^(r) plays 
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an essential role in the T-matrix element. Its calculation, how- 
ever, requires the wave functions of the two nuclei that are 
not readily available. We, therefore, substitute it by a con- 

1 /2 

stant & C xp times a physically plausible normalized function 
/(°) (r). In this work we take it to be a solution of 

(T r + y 17 (r)-e 17 )/(°)(r) =0, (11) 

where Vi 7 (r) is a Woods-Saxon potential that is supposed to 
represent the single-particle potential of p in 8 B. Since the 
norm of t^ ' (r) is the spectroscopic factor of p in 8 B, 6 eX p is 
its approximation if f^°'(r) is close to (^( ^(r)/||( / j(°^(r)|| 1 / 2 . 
After the substitution of (r), the T-matrix element be- 
comes 

<T ~ (2=l/2rr3-body 
T 3-body = ( XlX7 |C/ A7 + [/ A1 +y i7 |vI/ 3 . body ), (12) 

where 

* 3 -bod y = ^ (+) / (0) (r)e lP - R (13) 

is the wave function of the total system with incident wave of 
8 B consisting of a proton and an inert 7 Be bound in a state 
f( ' (r) impinging on the target with momentum P. 

Thus, in the particular case under consideration, the exact 
T-matrix element can approximately be calculated by three- 
body CDCC. Furthermore, since the reaction is peripheral in 
r as already mentioned, the radial part of f(°\r) in <J 3 ~ bod y 
can be replaced by its asymptotic form 

/(°)(r) ~ aW_^ 3/2 (2k r)/r. (14) 

Thus, 1 is approximately proportional to 

C = el£a. (15) 

1 12 

From the definitions of S CX p and a, it is clear that C is in 
good approximation the ratio of (pi°' (r) to W_ v 3 / 2 (2kQr)/r, 

i.e. the ANC of (p^ (r) we look for. Thus, one sees that 
the cross section is indeed proportional to C 2 . It should be 
noted that this proportionality holds only when the reaction is 
peripheral in r, which is indeed the case for the 8 B breakup 
concerned. 

In practice, we calculate the cross section (7 C aic from 
ij3-body f or eacn cjjojcg G f j(o) (y) w j m me three-body model 

CDCC as in the previous work 1 15, We then obtain 6 CXD 
by taking the ratio of the measured cross section to <j ca \ c . 
Since the amplitude a is readily known, the ANC, C, can 
then be calculated by Eq. dl5> . It should be noted that, since 
j(°)(r) and <p(°> (r) may have quite different forms in the in- 

1 /2 

ner region, 6 CX p may be quite different from the true spectro- 
scopic amplitude, H^ ^!")!! 1 / 2 , and depend strongly on the 
choice of /°(r). In contrast, C given by Eq. d!5l > is inde- 
pendent of / (0) (r) and exactly the ANC of </? (0) (r) if the 8 B 
breakup reaction is genuinely peripheral with no contribution 
of the inner part of the 8 B wave function. Once C is well de- 
termined from a peripheral reaction, Si 7 (0) can be obtained 
by 

5i 7 (0) w 38.0(1 - 0.0013a s )C 2 , (16) 



where a s is the s-wave p-'Be scattering length that can be cal- 
culated by the p- 7 Be interaction potential used in the analysis. 
The accuracy of the Si 7 (0) extracted with the ANC method 
can quantitatively evaluated by the dependence of C on the 
choice of (r), as shown in Sec. lIIICl 

Thus, we see that <Si 7 (0) can be extracted from the CDCC 
analysis of 8 B breakup based on the p+ 7 Be+A three-body 
model, if one concentrates on the 8 B breakup into forward 
angles, i.e. 8 S < 4°. 



B. E-CDCC and hybrid scattering amplitude 

In this subsection we describe the formalism of E-CDCC 
and how to construct a hybrid scattering amplitude from the 
results of E-CDCC and fully quantum-mechanical CDCC. 
The following formulation contains two important extensions 
of E-CDCC. One is the inclusion of the intrinsic spins of 
the projectile and its constituent fragments and the other is 
a more accurate treatment of the Coulomb wave function in 
the coupled-channel equations than in the previous formula- 
tion of E-CDCC jH. Detailed formalism of CDCC and its 
theoretical foundations are given elsewhere 

As already mentioned in Sec. Ill Al the total wave func- 
tion "fa-body of the p+ 7 Be+A three-body system satisfies the 
three-body Schrodinger equation, 

(iTi-body - T in ~ ei 7 )^ 3 -body(R, f) = 0. (17) 

In this model, the 8 B nucleus is assumed to consist of a 7 Be 
core in its ground state and a proton moving around it. Its 
intrinsic wave function, therefore, has the form $i^sim( r ) 
where £ is the orbital angular momentum of relative motion 
between p and 7 Be, S is the channel spin, and / and m are, 
respectively, the total angular momentum of 8 B and its z- 
component. The z-axis is taken to be parallel to the inci- 
dent beam. The index i stands for the i th of the states with 
the quantum numbers {£, S, I, to}. For simplicity we de- 
note below the five channel-indices {i, £, S, I, m} as c. The 
set {^c} consists of bound-states and "discretized-continuum- 
states" as in Refs. I15j [l a|. F or the calculation of the latter, the 
average metho d I l5lll6lll8l l36ll and the pseudo-state meth- 

^^^^1 ^^^3 ^^^3 ^^^^T^^ 

ods 1 151|37l 123139] are widely used. In the present work, we 
assume that the <f> c satisfy 

(Mr)|fc|« c (r)), = cmsi^c, (18) 

where h is the internal Hamiltonian of 8 B given by h = T r + 
Vi 7 (r), and ti^si is the intrinsic energy corresponding to the 
channel {%, £, S, I}. We take the following scheme of angular 
momentum coupling for < i ) c (r) 

$ c (r) = <fi,£si{r) ^2 (^meSms\Im)i e Ye me (r)Cs ms 

mi ,ms 

= Vi,isi(r) \i e Y t {T)®Cs] Im , (19) 

where Yi mi is the spherical harmonics and £sms i s tne s Pi n 
wave function of 8 B with channel spin S. 
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We expand ^ 



3-body 



in terms of the $ c (r) as 



* 3 -body(R, r) = J2 < fc(r)e-^ m " m °)^xc(i?, Or), (20) 



where the coefficient Xc(R> Or) exp[— i(m — mo)(j)ft\ of the 
expansion describes the cm. motion of the projectile rela- 
tive to A in channel c, and mo is m in the initial state. Since 
the $ c are chosen to form an approximate complete set in a 
finite region of space that is important for the reaction con- 
cerned [37], the expansion of Eq. i2Q\ is accurate in that re- 
gion of space. Unknown coefficients \c of the expansion are 
obtained by solving coupled-channel equations derived below. 

Multiplying Eq. (I17> by $*/(r) from the left, integrating 
over r, and summing over the spin states, one obtains the cou- 
pled equations for the Xc 

e- im '**(T R + a^esT - E in - e 17 ) Xc >(R,0 R ) 



l Xc(R,0 R ), (21) 



where Tr is the reduced kinetic energy operator defined by 



1 d / d\ d 



bdb \db) ' dz 2 



h (m' — mo)' 
27* b 2 



(22) 



in the cylindrical coordinates R = (z,b,(f> R ). The impact 
parameter of the collision, b, is defined by b — x 2 + y 2 
with R = (x, y, z) in the Cartesian coordinate. The two ar- 
guments of Xc in Eq. \2\\ are given by R = V z 2 + b 2 and 
Or = cos - 1 (zj \J z 2 + b 2 ) . In Eq. (I22t (j, is the reduced mass 
of the 8 B-A system. In the peripheral region, in which b is 
large, we neglect the last term of the r.h.s. of Eq. J22b . which 
turns out to have no effect on numerical results in the present 
case. Thus, 



Tr 



2[i 



-— (b— 

bdb V db 



dz- 



(23) 



The form factor F c i c in Eq. Mil of the coupling between chan- 
nels c and c' is given by 



F C , C (R) = <$ C '|C/A7 + C/Al|$ c ) r 

= F c , c (R,0R)e- i(m '- m) * R , 



(24) 



where we assume Ua7 and Uai to be central potentials. This 
assumption is valid as discussed in Ref. 1 30 ] since we use E- 
CDCC only to obtain the scattering amplitude for large values 
of b. The multipole expansion of T c i c is obtained from those 
of the coupling potentials 

U A7 (R A7 ) = ]>>i A 7 W/8)§£^(R)n,(r), 
U a1 (Rai) = J2 u AhR,rr/8)^J2 Y ^ Y ^)- 

A A v 



Then, T c i c can be expanded as 

T dc (R,Qft) = ^/- r V47(-) 2/ ' +m - s 



/(2l+l)(2l' + l)(2J+l)(2J' + l) 

x V W+W 

x (I -ml'm'\\m'-m){£0£'Q\\Q) 

x W{U'II'; \S)8 s >s 

x C\ (cos Oft) 

x / ri tl , s , r {r)UW{R,r)^si{r)r 2 dr, 



where W {it'll 1 ; XS) is the Racah coefficient, 



Cx ){ u + H)/2 / 2A + 1(A-H)! 

and (R, r) = u$ (R, r/8) + (R, 7r/8). 

Now we make the Coulomb-eikonal approximation 1 40] to 

Xc- 

X c(R,0ft)^i; c (b,z)cj>^+\b,z), 



(25) 



where (ft c is the Coulomb wave function with outgoing 
scattered wave and %jj c is the function to be determined with 
the Eikonal approximation. The Coulomb wave functions are 
given by 

<^c (+) (M = S^r(l+^ c )e^' R F£(R), 



(2tt) 3 / 2 ' 
G{-ir) c ,l,i(K c R-K c --R)) 



with r the Gamma function and G the confluent hypergeo- 
metric function. In actual calculation we make use of the ap- 
proximate asymptotic-form of 

1 



aC(+) 



C C(+) (M) 



(2tt) 3 / 2 



1 



i(K c R -K c R) 



xc 



(K c -R+?j c In [K e R-K c -H)) 



1 



(2tt) 3 / 2 



J(K c z+ Va In (K C R— K a z)) 



(26) 



which is valid for large values of b. The boundary condition 
for ip c in Eq. (125 \ is lim ip c (b, z) = <5 c o, where denotes the 

z — > — oo 

incident channel, so that Xc satisfies the appropriate boundary 
condition lim Xc(R, Or) = d ca ^ (+) (b, z). 

z — > — oo 

We make the following two approximations as in Ref. 1 30]. 
(i) For sufficiently large b, the flux of (f> c y ' is parallel to the 
z-axis. (ii) Local semiclassical approximation |41] to <fi c 



dz 



^Ki Kc (R)<f>^\R), 



(27) 



which is valid for a wave function in a slowly-varying poten- 
tial. The local momentum K C (R) is defined by 

^-Kl{R) = E ia + e 17 - e htS i - (28) 
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where Z%e (Zp^e) represents the charge of 8 B (A). Further 

approximation to cf> c was made in Ref. [30], i.e. the in- 
dependence of F£ was neglected. In the present paper, how- 
ever, we take account of that dependence using Eq. (I26> . This 
treatment of the Coulomb wave function slightly improves 
the agreement between the quantum-mechanical and eikonal 
scattering-amplitudes for large orbital angular momentum L 
of 8 B-A relative motion. The explicit forms of the two ampli- 
tudes are given later in Eqs. J37l >— J40i. 

With Eq. \23\ and the Coulomb-eikonal approximation, we 
obtain 



r RX J c {R,e R ) « t^(V c (M)0 c c(+) (M) 



ih 



+ (t^ c c (+>(M)W(M 



(29) 



where the second-order derivative of ip c is neglected since it 
is slowly varying compared with <pc . Inserting Eqs. MAI . 
(|26}, {Zfr, and 



Z s Z A e 2 



R 



into Eq. 12 H . we arrive at the following coupled-channel 
equations 



ih 2 



J(K C ,-K e )z 



(30) 



where the non-dynamical variable b is put in a superscript, the 
reduced coupling-potential $ cc > is defined by 



and 



irj c i In (K C ,R-K c ,z) 
piri c \n{K c R-K c z) 



R 



Sec', 



{K C ,R-K d z)^' 
(K C R - K c z) 



17),; 



Since Eq. d30i is a set of first-order differential equations and 
contains no coefficient with very large angular momentum L, 
one can solve it with high computational speed and accuracy. 

The eikonal scattering amplitude with Coulomb distortion 
is given by / C E = f? xlth 5 c0 + /c0 > where /c0 Uth is the Ruther- 
ford amplitude and 

(2tt)V 



JcO — 



H 2 



/£^(*)*f ( - } (M)*? +) M 



xe 



^V{z)dR 



(32) 



with 



<tP c t-\b,z) 



(2tt) 3 / 2 
1 

(2tt) 3 / 2 ' 



r^i + ^e^ 11 ^ (R) 



JK'c-U -irjc In (if jR+K^R)) 



(33) 



The outgoing cm. momentum J£' c is chosen to be on the z-x 
plane, following the Madison convention. For forward scatter- 
ing in which K.' c is almost parallel to K c s i.e. to the z-axis, we 
approximate, in the same way as in Ref. 1 30], the phase-factor 
(K c ; - K' c ) ■ R, which appears in < />* C( " ) (6, z)^, (+) (b, z), 
by — K c 9fb cos (f> R + (K c i — K c )z, where 9f is the scattering 
angle of 8 B, and K C R + K' c • R in Eq. 03 by K C R + K c z. 
Under these approximations, 



IE 
cO 



2irh 2 



-i(m-m )<t>R e -iK c sin f b cos tf> R 



x W W / ^ (z) ^ (z) e^' -"J"*® (z)dz 



(34) 

with H c b) = exp[2i?7 c ln {K c b)\. 

The integration over z in Eq. J34l > can be done analytically 
with the help of Eq. POL with the result 



J c< 



ih 2 



K c (R)4 b \z) 



-K r 



S. 



(b) 
rO 



(35) 



where dK c (R) / dz is neglected, since K C (R), with large b, 
is a slowly- varying function of z. The eikonal S'-matrix el- 
ements are defined by S^q — \\m%jj^\z). Insertion of 



(31) Eq. O into Eq. <03J leads to 



f/E ^ Kc 
/c0 ~ 2ni 



'Hi b) 



^ — i(m — 7no)<pR^ — iK c sin 9 fb cos (j>n 



0(6) i 



bdbd(f) R 
(36) 



We now transform the integration over b to the summation 
over L as in Ref. l3(ill . The resulting eikonal scattering ampli- 
tude is 



f E _ --Ruthr , ^ WE v 

JcO - JcO °cQ + -77- JL-cO Yl ' 



,(K'), (37) 



where 



JL: 



L:c0 



H (bo;L 



2L+1 
47T 



(m— mo) 



5. 



(38) 
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with K c b C ;L = L + 1/2. The quantum mechanical scattering 
amplitude obtained by CDCC is given by 



/, 



<::0 



Ruth c 
JcO °c0 



2tt x 



f' L Q c0 Y Lm „ mo (K') (39) 



with 

f 'Q 

J L:cO 



2L + 1 
47T 



L+7 J+/ 

E E 

J=|£-J| i =|J-Io| 

x (IomoLoO\ Jmo)(ImL mo—m\Jmo) 

X ( S iLtSI,i L i S I ~ S iio S LL Su 5sS 5lI ) 
xe i(n+^ )(_)m-ra 0] (4Q) 



where <tl is the Coulomb phase shift and J is the total angular 
momentum of the three-body system. We define the hybrid 
scattering amplitude by 

f H _ f Ruth c , 2?r \^ f /Q v /x>/\ 
JcO — JcO °c0 ' JL;cO ■ r im-m V rv J 



L<L C 



~j7~ E ^im-mofK') 
* ° L>L C 



(41) 



and use it in the analysis of 8 B breakup. 

The connecting-angular-momentum Lq between quantum 
mechanical and eikonal amplitudes is chosen so that ff^. cQ 

agrees with /^. c0 for L > Lq. It should be noted that Eq. (14- 1 i 
includes all quantum effects and also the interference between 
amplitudes in the two L-regions. This hybrid CDCC calcu- 
lation turns out to be capable of evaluating the nuclear and 
Coulomb breakup amplitudes with very high accuracy and 
computational speed. In fact, we have obtained the same ac- 
curacy as standard fully quantum mechanical CDCC calcula- 
tions in about 1/10 of computing time. 



III. RESULTS AND DISCUSSION 



A. Numerical calculation 

For the wave functions of 8 B that form the basis of the 
CDCC expansion of the three-body system we make the fol- 
lowing assumptions. We assume that the ground state wave 
function of 8 B, with the z-component too of the total spin 
Iq = 2 of it, is given by 

/W(r) - { Wl [»yi(r) ® Ci] 2mo + w 2 [^(r) ® C 2 ] 2m J 



x/(°>(r) 

^i/i ( l(r)+^2/ 2 ( l(r) 



(42) 



in the channel spin representation, neglecting the £ ^ 1 com- 
ponents. The values of the coefficients w\ ~ 0.397 and 
W2 ~ 0.918 are obtained from the result of the microscopic 



calculation |17J| of 8 B, Cf /2 /C. 



y 3/2 



0.157, where C 1/2 



(C3/2) is the ANC in the proton spin-orbit coupling represen- 
tation corresponding to the j = \l+ ^\ = \/2 (3/2) state. We 
further assume the effective p- 7 Be interaction for the ground 
state of 8 B, to be independent of S, hence the same ra- 
dial dependence (r) of the S = 1 and S = 2 components 
in Eq. J42I . This assumption is valid since only the tail region, 
r > i?N, of f(°^ (r) contributes to the 8 B breakup process con- 
cerned. The wave function of the p- 7 Be scattering state with 
outgoing scattered waves is given by 



/£l(k,r) 



is 



£17 -T r - Vvi 



l£ 



r Cs, ms , 5=lor2, 
(43) 

where e lk r Cs,m s i s trie incident wave with definite channel 
spin S and its z-component mg, £17 is the incident energy in 
the cm. frame of the p- 7 Be system, and k is the incident mo- 
mentum. The wave function with incoming scattered waves 
fs~m (kj r ) corresponding to /g^ (k, r) used below is de- 
fined by (k, r) = f <% g (-k, r). 

We assume that V17 has no tensor force component. We 
consider the breakup of 8 B to its scattering states at low £17, 
and truncate the modelspace by £17 < 10 MeV. In the exper- 
iment we analyze, however, the cross section is measured for 
£17 in the range of 500 keV < £17 < 750 keV, and it turns out 
that only coupling of those states with £17 much lower than 
lOMeV need be taken into account as described in Sec. lIIlBl 
for which only s-waves are affected by the nuclear part of V17. 
All the higher partial waves have negligible amplitude inside 
its range because of the Coulomb and centrifugal barriers and 
there have no coupling with the s-waves. Therefore, all the an- 
gular momenta £, S and / are conserved in /^^(k, r). The 
channel spin S is conserved also in s-waves since S — I that 
is a good quantum number. Thus, one sees that fg^ (k, r) is 
the pure state of S, while £ and / are mixed as in the incident 
wave e jk r . It should be noted, however, that Vn for the s- 
waves depends on S as is evident from the large S'-dependence 
of the s-wave p- 7 Be scattering lengths: af =2 = —7 ± 3 fm 
anda^ =1 = 25±9fmEl. 



Furthermore, it turns out that the spin-dependent part of the 
p- 208 Pb optical potential has no effect on the resulting elastic 
cross section at forward angles. Thus, the channel spin S is 
conserved during the breakup process of 8 B by the 208 Pb tar- 
get fiUl . Under these circumstances, the S'-dependent 3-body 
T-matrix element in actual CDCC calculation, ^g^^, is 
given by 



-3-body 



(-) 
Sms 



ws(f, 
x|fi(+>/, (0) 



'k, r)e 



PR 



U A 7 + U A 1 



S, mo 
,,, tr-'3-body 



(r)e 



PR\ 



(44) 



where is the wave matrix given in Eq. 

The unpolarized 8 B breakup cross section calculated by 
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CDCC, a, is given by 
Vp 



a= 2I + 1 E \^2{h m ^l m r\Sms)ws%l 



'3-body 

msma 



2/n 



momim 7 Sm s 

"S,rn 5777,0 



^ „, 2 l^'3-body 2 



T E 



E 



(45) 



where mi (7717) is the z-component of the spin of p ( 7 Be) 
and the orthonormality of the Clebsch-Gordan coefficient is 
used. We compare a with the unpolarized cross sections mea- 
sured by RIKEN experiment |7J, |8j, |9| . In the calculation of 
<7i and (T2, we use an /-independent p- 7 Be effective inter- 
action Vn as mentioned above. Then, the states of 8 B with 
fixed 5 and £ and different values of I are degenerate, which 
greatly simplifies numerical calculation. The value of a s in 
Eq. (116 ) corresponding to the spin-dependent ANC analysis 



is Il34ll a s = u> 2 af _1 

r(0) 



„2„S=2 



-2.8 fm. 



For V£>, the single-particle potential of Esbensen and 
Bertsch (EB) 1E3I1 is used, except that we neglect the spin- 
orbit part of the potential and adjust the depth of the central 
potential to reproduce the separation energy of p, 137 keV. 
The scattering states for the p-waves are also calculated with 
this vffl . For d- and f-states with both 5 = 1 and 5 = 2, we 
use the potential of Barker 14411 . For the s-state in the 5 = 2 
channel, we use the potential of Barker that gives the s-wave 
p- 7 Be scattering length af =2 of —8 fm. For the s-state with 
5 = 1, we change the depth of the potential of Barker to 25.7 
MeV so that the resulting scattering length af =1 is 25 fm. 

As for the distorting potential between p ( 7 Be) and 208 Pb 
we adopt the global optical potential by Koning and De- 
laroche 14511 (Cook 1 46]). We neglect the spin-orbit parts 
of the p- 208 Pb potential. The multipoles for nuclear and 
Coulomb coupling-potentials are included up to A = 6. The 
discretized-continuum states of 8 B are constructed by the av- 
erage method 1 15J [Kj []Jl l3al . The maximum excitation en- 
ergy of 8 B is 10 MeV and 10, 20, 10, and 5 discretized- 
continuum states are taken for the I = 1, 0, 2, and 3 states, 
respectively. The resulting number of scattering-channels is 
138. The maximum values of r, R, and L are, respectively, 
200 fm, 1000 fm, and 12000. The modelspace described 
above turns out to give good convergence of the resulting 
breakup cross sections. 



B. Analysis of B breakup experiment 

Figure 2 shows the convergence of the hybrid calculation 
with CDCC and E-CDCC for the breakup cross section of 8 B 
by 208 Pb at 52 MeV/nucleon. The cross section integrated 
over the excitation energy of 8 B measured from the p+ 7 Be 
threshold energy, £17, from 500 keV to 750 keV is plotted 
as a function of 9$. The dashed line is the result of the E- 
CDCC calculation (Lq = 0) and the dotted, dash-dotted, and 
solid lines correspond to the hybrid calculation with Lq — 
500, 1000, and 1500, respectively. The figure shows that the 
hybrid calculation converges with Lc = 1000. The difference 
between the dash-dotted and solid lines is only less than about 



1% in magnitude. Thus, we regard the hybrid calculation of 
CDCC and E-CDCC with L c = 1000 as the fully quantum- 
mechanical CDCC calculation. 

One sees in Fig. 2 the oscillation of the cross section at 
forward angles around 2°, which turns out to be due to the in- 
terference between nuclear and Coulomb breakup-amplitudes. 
Therefore, it can be concluded that nuclear interactions affect 
the breakup cross section even at very forward angles. How- 
ever, as one sees below, this oscillation is not observed in ac- 
tual experimental data because of the limit of the resolution of 
6>8 in the measurement. 

Figure 3 shows the comparison between CDCC calcula- 
tion and the experimental data for the 8 B breakup. In order 
to take account of the experimental resolution and efficiency, 
the theoretical result has been smeared out by using the filter- 
ing table B47I1 provided by the authors of Ref. |8]. The upper 
panel shows the filtering effect on the cross section calculated 
with CDCC, where the solid and dashed lines represent the re- 
sults obtained with and without the smearing procedure. In the 
lower panel, the smeared total breakup cross section is decom- 
posed into the s- (dashed line), p- (dotted line), d- (dash-dotted 
line), and f-state (dash-two-dotted line) breakup components 
of 8 B. 

One sees the smeared result of CDCC well agrees with the 
experimental data at forward angles (0g < 4°), while it under- 
estimates the data at backward angles. The lower panel shows 
the importance of the p-state breakup for 9% > 5°, which im- 
plies that nuclear breakup is significant in this angular region. 
Thus, careful description of the p-state 8 B wave function, the 
resonance structures of it in particular, will be necessary to re- 
produce the experimental data at backward angles. The under- 
estimation there may be due also to dynamical roles played by 
the excited-core component in 8 B, 7 Be(l/2~)£g)p(3/2 _ ), that 
is neglected in the present work as discussed in Sec. Ill Al Al- 
though further study in this line will be very interesting, it is 



Eikonal(L c = 0) 

Hybrid (L c = 500) 

Hybrid (L c = 1000) 

Hybrid (L c = 1500) 




500keV<s 17 <750keV 



2 4 6 

6 8 (deg) 



8 10 



FIG. 2: The B breakup cross section integrated over £17 in the range 
of 500-750 keV, as a function of #8- The dotted, dash-dotted, and 
solid lines represent the results of the hybrid calculation of CDCC 
and E-CDCC with L c = 500, 1000, and 1500, respectively. The 
result of E-CDCC, i.e. Lc = 0, is also shown by the dashed line. 
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beyond the scope of the present paper. Additionally, the fil- 
tering table used was made with assuming the s-state breakup 
of 8 B. Thus, quantitative comparison between the calculation 
and the experimental data, to extract 5*17(0) with high accu- 
racy, can only be done in the region where the s-state breakup 
cross section is dominant. Therefore, below we use the data 
for 6g < 4° to determine 6 oxp by 



where a is given by Eq. ( 145 1 and cr cxp is the corresponding 
experimental data. 

Figure 4 shows the result of the x 2 -fit of the theoretical cal- 
culation to the experimental data. The solid line is the same as 
that in each panel in Fig. 3 but multiplied by the spectroscopic 
factor 6 oxp = 1.11. The dashed line represents the result 
before being fitted, i.e. with S cxp = 1.00. Each horizontal 
bar put on the eight data points below 4° does not represent 
a statistical error but it shows the range of 6*§ in which the 
breakup cross sections contribute to each data point The 




2 4 6 8 10 




6 8 (deg) 

FIG. 3: (a) The calculated 8 B breakup cross sections with (solid line) 
and without (dashed line) the resolution and efficiency of the mea- 
surement [47] taken into account. The experimental data are taken 
from Ref. (b) Decomposition of the smeared 8 B breakup cross 
section into the s- (dashed line), p- (dotted line), d- (dash-dotted line), 
and f-state (dash-two-dotted line) breakup components. 



1000 




6 8 (deg) 

FIG. 4: The result of the x 2 - nt of the breakup cross section calcu- 
lated with the 8 B single-particle potential of Esbensen and Bertsch 
(EB) l43ll . The solid line is the same as that in each panel in Fig. 3 
but multiplied by the spectroscopic factor 6 cxp = 1.11. The result 
with S cxp = 1.00 is also shown by the dashed line for comparison. 
For the eight data points below 4°, which are used in the fitting pro- 
cedure, we have put the horizontal bar that represents the resolution 
of 8s 1 8] of the measurement. 



value of x 2 P er datum obtained is 8.5. Although the qual- 
ity of the fit is not good, it should be noted that estimation 
of error-propagation in the present case, i.e. comparison be- 
tween smeared experimental data and a smeared numerical 
result, is very complicated; the value of \ 2 per datum shown 
above takes account of no error with respect to 6* 8 . The value 
of a in Eq. (fT4b for the EB 8 B wave-function is 0.704 fmT 1 / 2 . 
With the values of 6 oxp and a, the ANC C is obtained from 
Eq. <BJ as C = 0.740 fm" 1 / 2 . One can then determine 
Si 7 (0) by inserting this result into Eq. dl6> . together with 
a s = -2.8 fm, i.e. Si 7 (0) = 20.9 eV b. In the next subsec- 
tion we quantitatively evaluate uncertainties of the extracted 
Si 7 (0). 



C. Uncertainties of the extracted Si7(0) 

First, we evaluate the uncertainty of Si7(0) that comes 
from the use of the ANC method, just in the same way as in 
Refs. We show in Fig. 5 the result of CDCC calcula- 

tion with the p-state single-particle potential of 8 B by Kim et 
al. 1 48 ] (the dashed line) and that multiplied by S cxp = 0.867 
(the solid line); the latter agrees with the result shown in Fig. 4 
within 1%. One sees that S cxp indeed depends on the p-state 
potential. Actually, a in Eq. dl4t for Kim's potential is 0.796 
fm -1 / 2 , which is quite different from the value for the EB 
potential, 0.704 fm -1 / 2 , because of the difference in the ge- 
ometry of the two potentials of about 20%. On the contrary, 
the value of C in Eq. d!5l > calculated with Kim's potential is 
0.741 fm" 1 / 2 that agrees very well with the result obtained 
with the EB potential shown in Sec. lIIlBl This result shows 
that the ANC method works with very high accuracy in the 
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FIG. 5: Same as in Fig. 4 but for the CDCC calculation with Kim's 
single-particle potential 1 48] for the p-state of 8 B; the spectroscopic 
factor Scxp obtained is 0.867. The result with Kim's potential before 
X 2 -fitting, namely with S C x P = 1.00, is also shown by the dashed 
line for comparison. 



FIG. 6: Dependence of the B breakup cross section on the s- 
state single-particle potentials of 8 B for the S — 1 and S = 
2 channels. The solid line represents the result with the choice 
of the potentials that corresponds to the s-wave p- 7 Be scattering 
lengths {af =1 ,af= 2 } = {25 fm, -8fm}. The dashed (dotted) 
line shows the result with the s-state potentials corresponding to 
{of =1 , af= 2 } = {16 fm, -10 fm} ({34 fm, -4 fm}). 



present analysis, i.e. the error of the ANC method is negligi- 
ble. 

Second, we estimate the effect of ambiguity of the distort- 
ing potentials used in the CDCC calculation on the obtained 
Si7(0). For this purpose we use alternative optical poten- 
tials forp- 208 Pb and 7 Be- 208 Pb. For the p- 208 Pb potential, we 
modify the parameter set of Koning and Delaroche so that the 
calculated elastic cross section without spin-orbit terms repro- 
duces the one obtained with the full components of the optical 
potential. The parameter set thus obtained is Vy = 39.82 
MeV, W v = 0.776 Me V, and W D = 13.47 MeV, with the 
same notation as in Ref. [45]; all other parameters are not 
changed. For the 7 Be- 208 Pb potential, we use a single-folding 
model with the nucleon- 208 Pb potential of Koning and De- 
laroche, neglecting the spin-orbit terms. The density distribu- 
tion of 7 Be is assumed to be of the Gaussian form with the 
ran ge t hat reproduces the rms matter radius of 7 Be of 2.48 
fm I49B . The difference between the results calculated with 
the original and alternative optical potentials is found to be 
negligibly small (not shown). We conclude, therefore, that the 
error of £17(0) that comes from the ambiguity of the distort- 
ing potentials in the present case is negligibly small. 

Finally, we estimate the uncertainty of 617(0) due to the 
ambiguity of the s-state single-particle potentials of 8 B for the 
5 = 1 and 5 = 2 channels by changing the depth of the 
potential for the S = 1 (S = 2) channel ranging from 23.3 
MeV (52.8 MeV) to 30.5 MeV (57.6 MeV), which gives the 
range of the resulting s-wave p- 7 Be scattering length af =1 
(af =2 ) from 16 fm (-10 fm) to 34 fm (-4 fm); these values 
of of =1 and of =2 are within the range of errors of the ex- 
perimental values of them l42Tl . In Fig. 6 the solid, dashed, 
and dotted lines represent the results with the s-state po- 
tentials corresponding to {af =1 , af =2 } = {25 fm, — 8 fm}, 
{16 fm, — 10 fm}, and {34 fm, — 4 fm}, respectively. One 
sees the breakup cross section slightly depends on the choice 



of the s-state potentials of 8 B; the resulting £17(0) is found to 
vary from 20.3 eV b to 21.9 eV b. Thus, we evaluate the error 
of Si 7 (0) to be +4.9%/-2.9%. 

Summarizing the discussion given above, the result of the 
present paper is as follows. The central value of £17(0) is 20.9 
eV b. The theoretical error of the extracted £17(0) concerned 
with the s-wave p- 7 Be scattering length is +4.9%/— 2.9%. Af- 
ter including the 8.4% systematic experimental error, we ob- 
tain 

£17(0) = 20.9t£;g (theo) ± 1.8 (expt) eVb. 



D. Discussion of the extracted 517(0) 

The main result of the present paper is 5i 7 (Q) = 20.9l Q 'g 
(theo) ±1.8 (expt) eV b derived from 208 Pb( 8 B, p+ 7 Be) 208 Pb 
at 52 MeV/nucleon. This value is significantly larger than 
5 17 (0) = 18.9 ± 1.8 eV b [9] obtained in the previous anal- 
ysis of the same experiment with the first-order perturbation 
theory. In order to clarify the reason for the difference, we 
discuss roles of nuclear interaction, E2 transitions, and higher- 
order processes. In Fig. 7 the solid line is the same as in the 
upper panel in Fig. 3 that shows the result of full CDCC. The 
dashed line corresponds to the calculation without breakup 
components of nuclear coupling potentials. It should be noted 
that their diagonal components are included in the calculation 
to take account of the absorption of the flux of the incident par- 
ticle by the target nucleus. One sees that the dashed line agrees 
with the solid line below 2° but deviates from it for 9s > 2°. 
As a result, the value of Si 7 (0) obtained, Si 7 (0) = 20.0t} g 
eVb, is smaller than 5i 7 (0) = 20.9±J;g (theo) ±1.8 (expt) eV 
b obtained with the full CDCC calculation mentioned above. 
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FIG. 7: Same as in Fig. 3 with different assumptions for CDCC. 
The lines represent the results of CDCC with nuclear and Coulomb 
couplings and all higher-order processes (solid line), CDCC without 
nuclear breakup (dashed line), one-step CDCC with Coulomb El and 
E2 components (dotted-line), and one-step CDCC with Coulomb El 
only that corresponds to the first-order virtual photon absorption the- 
ory (dash-dotted line). 



The dash-dotted line in Fig. 7 represents the result of first- 
order iterative solutions of CDCC, designated as one-step 
CDCC, without nuclear breakup, including only A = 1 (El) 
component of the Coulomb interaction; this calculation is es- 
sentially the same as the first-order perturbation theory (vir- 
tual photon theory) used in the previous analysis of the exper- 
imental data 0. It overestimates the solid line above 1.5° and 
the resulting value of 5*17(0) is 19. 0~t\'^ eV b, which agrees 
well with the value 18.9±1.8eVb obtained in Ref. H- 

If one includes the A = 2 (E2) component in the one- 
step CDCC calculation of Coulomb breakup, the dotted line 
in Fig. 7 is obtained; the coupling potentials with A = 2 
are artificially multiplied by 0.7, following the analysis of 
the MSU data 111 311 . One sees from Fig. 7 that inclusion of 
the reduced E2 component somewhat increases the breakup 
cross section, which results in further decrease of 5i7(0) to 
Sn(0) ~ 17.9 jl^'g eV b. This result is consistent with the 
conclusion of Ref. Q, in which 5i 7 (0) = 17.8±\-j eV 
b was derived with first-order perturbation theory including 
both El and the reduced E2 components. If the E2 compo- 
nent is not scaled, the resulting value of £17(0) is 16.7 eV b. 
This value is about 20% less than the result of full CDCC, 
20.9 eV b. The difference is due to the nuclear and higher- 
order Coulomb processes since the coupling potentials with 
A > 3 are found to have little effect on the total breakup- 
cross-section at #g < 4°. 

Thus, description of 8 B breakup process with nuclear and 
Coulomb breakup processes of both the El and E2 transitions 
and higher-order processes is a key to solve a puzzle recog- 
nized so far of the discrepancy between the values of 5*17(0) 
extracted from direct p-capture reactions, 5i7(0) = 22.1±0.6 
(expt) ±0.6 (theo) eV b, and indirect 8 B dissociation experi- 
ments. In order to clarify the role of these components in 



8 B breakup reaction in general, it is necessary to carry out 
analyses to wider range of data, such as those measured at 
MSU EQ and GSI E3IH, including the data on other 
quantities than that dealt with in the present work. Analysis 
of parallel momentum {pi 1) distribution of 7 Be-fragment after 
breakup of 8 B is particularly important, since the role of E2 
component was determined from it II 2ul 311. In the analysis in 
Ref. lll3ll the result of first-order perturbation theory with an 
adjustable parameter was compared with CDCC calculation 
that is essentially the same as in the present paper. Because of 
the agreement of the two results, it was concluded that higher 
order processes were unimportant without quantitative evalu- 
ation of the nuclear and higher-order Coulomb processes as 
done in the present work. Actually, the contribution is about 
20% of the total as already described. A CDCC analysis of 
the pn distribution has been done by Mortimer et al. l28ll . 
where the dependence of the calculated result on the strength 
of the Coulomb A = 2 (E2) coupling potentials was discussed. 
They found that the E2 potentials must be multiplied by 1 .6 to 
reproduce the experimental data, although the origin of the 
enhancement was not clear. Therefore, further investigation 
of the accuracy of the enhancement factor 1 .6 seems neces- 
sary. Unfortunately, filtering table for the MSU experiment 
is not available to the present authors. At this stage, there- 
fore, we can quantitatively extract 5i7(0) with CDCC and the 
ANC method only from the 8 B breakup experiment done at 
RIKEN. Nevertheless, comparison of our pure theoretical re- 
sult with that shown in Ref. 1 28 ] is very interesting, which will 
be a subject of our future work. 



E. ANC analysis of the direct measurement 

In this subsection we report on our application of the ANC 
method to the p-capture reaction, 7 Be(p, 7) 8 B, at energies of 
1 16-362 keV |4]. We assume direct capture of p by 7 Be to the 
ground state of 8 B by a pure El transition. A simple potential- 
model is assumed for the bound and the scattering states of the 
p+ 7 Be system, using the single-particle potential of Sec. lIII Al 
The s- and d-state components in the initial scattering-state 
wave-function are taken into account, and the intrinsic spins 
of p and 7 Be are included in the channel spin representation 
described in Sec. IHI Al In Fig. 8 the solid line shows the cal- 
culated 5i7 as a function of the incident energy E c . m . of the 
reaction in the cm. frame of the p- 7 Be system . The plotted 
5i7(-E c .m.) has been multiplied by the spectroscopic factor 
S 0X p = 1.15 obtained by the x 2 -fit to the data below 362 keV. 
The calculated Sn(E c . nl .) is not reliable at high energies, say 
Ec.m. > 400 keV, since the inner part of the 8 B wave function 
plays an important role. At low energies, however, the ANC 
method works very well and a reliable value of the ANC C, 
and consequently 5i7(0), can be derived from the x 2 -fitting 
procedure. The value of 5i7(0) thus obtained is 21.7 eV b. 
Uncertainties of the 5i7(0) are evaluated in the same way as 
in Sec. lmCl and found to be +0.37/-0.24 eV b, almost all 
of which is due to the ambiguity of the s-state potentials for 
the 5 = 1 and 5 = 2 channels of the p- 7 Be system. In- 
cluding the 2.3% systematic experimental error, we obtain the 
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FIG. 8: The astrophysical factor Sn as a function of the incident 
energy -E c .m. of the reaction in the cm. frame of the p- 7 Be system. 
The solid line is the result of the calculation with a simple potential 
model, multiplied by <3 cxp = 1.15. The experimental data are taken 
from Ref. 0|. 



following result: 5i 7 (0) = 21.71^ (theo) ±0.50 (expt) eV 
b. This result is consistent with both the value extracted from 
8 B dissociation, S 17 (Q) = 20.9±J;g (theo) ±1.8 (expt) eV b, 
described in the previous subsections, and the result obtained 
in Ref. 0], 5 17 (Q) = 22.1 ± 0.6 (expt) ±0.6 (theo) eV b, al- 
though no evaluation of the ambiguity due to the uncertainty 
of the s-wave p- 7 Be scattering lengths, as the one described in 
Sec. Him is made in Ref. H . 



IV. SUMMARY 



partial waves. Improvements of E-CDCC are made to cope 
with the Coulomb distortions of the incident and scattered 
waves of relative motion and intrinsic spins of the projectile 
and the ejectiles. 

The main source of the error in the obtained 5x7(0) is the 
ambiguity of the s-wave p- 7 Be interaction potentials for the 
channel spin 5 = 1 and 5 = 2 states. The errors asso- 
ciated with the use of the ANC method and the ambiguity 
of the distorting potentials are both negligible. Calculations 

28 r 1 

ANC: 

26 - 21.7tS 55 eVb 

' \ FullCDCC: 

24 ~ T Junghans[4] \ 20.9^;" eVb 

i'T\ V\ One-step El: 

«, 22 r RIKEN19] i ,T i ; lon+1 o 



16 rMSU[13] A direct capture One-step E1-E2: 
[ • 8 B breakup 17.9!};^ eVb 

FIG. 9: The values of Si 7 (0) extracted by the ANC method are 
shown by open symbols: the circles are obtained from three-types of 
CDCC analysis of 8 B breakup and the triangle is the result of analysis 
of 7 Be(p, 7) 8 B with a simple potential model. The error of each 
result is obtained by adding the theoretical and experimental errors 
in quadrature. The results are compared with Sn(0) extracted from 
the precise measurement of direct capture (p, 7) cross section 0] 
(closed triangle) and those obtained from 8 B dissociation with first- 
order perturbation theory (closed circles). 
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eVb 
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19.0 



A direct capture (this work) 



o 8 B breakup (this work) 
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The principal result of the present paper is the value of 
5i 7 (0) of 20.9±J;g (theo) ±1.8 (expt) eV b obtained by 
an analysis of the cross section of the 8 B breakup reac- 
tion 208 Pb( 8 B j?+ 7 Be) 208 Pb at 52 MeV/nucleon measured 
at RIKEN (21 |H Lai by means of the method of continuum- 
discretized coupled-channels (CDCC) combined with 

the asymptotic normalization coefficient (ANC) method 1 17]. 
The value is consistent with the one extracted from the 
precise measurement of the cross section of direct capture 
7 Be(p, 7) 8 B, 5x7(0) = 22.1±0.6 (expt) ±0.6 (theo) eV b |4]. 

The CDCC calculation is based on the three-body model 
of p, 7 Be, and the target 208 Pb nucleus, with 7 Be staying in 
the ground state throughout the reaction process. This model 
is shown to be adequate in the breakup reaction in which the 
fragments p and 7 Be are ejected in very forward angles, less 
than 4°, which is the case for the experimental data analyzed 
in this work. Nuclear and Coulomb El and E2 transitions 
and multi-step processes of all-order are included in the cal- 
culation. The calculated cross section is shown to be propor- 
tional to the squared asymptotic amplitude of the overlap of 
the ground state wave functions of 7 Be and 8 B, which is a nec- 
essary condition for the applicability of the ANC method. The 
CDCC calculation is made efficient by the use of the eikonal 
CDCC method US (E-CDCC) for large angular momentum 



of 5i7(0) with some simplified assumptions are summarized 
in Fig. 9. The results of first-order iterative solutions of 
CDCC, designated as one-step CDCC, correspond to and 
agree with those of first-order perturbation theory in the pre- 
vious work 0, U [3 fl3ll • The inclusion of the Coulomb 
quadrupole (E2) transitions, scaled by 0.7, in the one-step cal- 
culation decreases 5i7(0) by about 6%, and multistep pro- 
cesses increase it by about 20%. This shows the crucial im- 
portance of accurate description of the 8 B breakup process by 
CDCC including nuclear and Coulomb El and E2 transitions 
and all higher-order processes. It will be interesting to apply 
the method of the present paper to analyses of wider range of 
experiments such as those in Refs. llOt fill fl2l fl3ll . It may, 
however, be that the three-body CDCC used in the present pa- 
per is not valid in general and the use of CDCC with four-body 
model such as the one in Ref. 1 50 ] becomes necessary. 
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